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Abstract. In this talk I will present a model for primordial galaxy for¬ 
mation. In particular, I will review the feedback effects that regulate the 
process: (i) radiative (i.e. ionizing and H 2 -photodissociating photons) 
and (ii) stellar (i.e. SN explosions) feedback produced by massive stars. 
I will show the results of a model for galaxy formation and IGM reion¬ 
ization, which includes a self-consistent treatment of the above feedback 
effects. Finally, I will describe a Monte Carlo method for the radiative 
transfer of ionizing photons through the IGM and discuss its application 
to the IGM reionization problem. 


1. Introduction 

The application of the Gunn-Peterson test to QSOs absorption spectra suggests 
that the Intergalactic Medium (IGM) is completely reionized by z ~ 6. Several 
authors claim that the known population of quasars and galaxies provides ~ 10 
times fewer ionizing photons than are necessary to keep the observed IGM ion¬ 
ization level; thus, additional sources of ionizing photons are required at high 
redshift, the most promising being early galaxies and quasars. Recent observa¬ 
tional evidences suggest the existence of an early population of pregalactic stellar 
objects which could have contributed to the reionization and metal enrichment 
of the IGM. In the following section, I will review the processes that regulate 
the formation and the evolution of such galaxies. 


2. Primordial Galaxy Formation and Feedback Effects 

Once the gas has virialized into the potential well of preexisting dark matter 
haloes, additional cooling is required in order for the gas to further collapse and 
form stars. Objects with virial temperatures T^ir > 10^ K rely on hydrogen Lya 
line cooling for their formation, while the main coolant for objects with Tyir < 
10^ K {PopIII objects) is molecular hydrogen, a key species for galaxy formation 
in the early universe. H 2 is a very fragile molecule, as it is easily dissociated by 
photons with energy in the range 11.26-13.6 eV. Thus, the existence of an UV 
flux below the Lyman limit, capable of dissociating H 2 in collapsing structures 
can strongly influence the formation of PopIII objects. Giardi et al. (2000) have 
derived the minimum mass that a forming galaxy requires to be able to shield 
against an external incident dissociating flux: they found that the mass is in the 
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Redshift 

Figure 1. Evolution of HII filling factor as a function of redshift 
(solid curve). The contribution from objects with masses in different 
range is also shown: 5 — 10 x 10 ®Mq (dotted line), 1 — 5 x IO^Mq 
( short dashed), 5 — 10 x IO’^Mq (long dashed), 1 — 5 x 10 ®Mq (short 
dashed-dotted), 5 — 10 x 10®M© (long dashed-dotted), 1 — 5 x 10 ®Mq 
( short-long dashed). 


range 10^Mq < M < 10®Mq, depending on the intensity of the incident flux 
and the formation redshift. 

Once an object has collapsed, the first episode of star formation can strongly 
influence the subsequent star formation process via the effects of mass and energy 
deposition due to massive stars. These processes may induce different phenom¬ 
ena, such as a heating and a loss of the interstellar gas, or a dissociation of H 2 
molecules in star forming regions. As a consequence, the star formation rate 
decreases, more or less dramatically depending on the mass of the object. 


3. IGM Reionization 

Ciardi et al. (2000) have studied the interplay between the primordial galaxy 
formation process and the IGM reionization. In their analysis the authors have 
included a self-consistent treatment of the feedback effects that regulate the 
formation and the evolution of galaxies, as explained in the previous section. 
One of the output of their model is the redshift evolution of the HII filling factor, 
full (Fig. 1). At 2; > 15 the main contribution to fnii comes from small mass 
objects (M ~ IO^Mq), while at lower redshift, when the formation and evolution 
of such objects is suppressed by feedback effects, the main contribution comes 
from galaxies of M > 10 ®Mq. Thus, the contribution of small mass objects 
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Figure 2. Ionized regions in a 30 Mpc comoving simulation box at 
z ~ 13 (left panel) and z ~ 10 (right panel). 


to the ionizing photon budget is negligible. These galaxies are nevertheless 
important as sinks of radiation on small scales. 

To improve on the mentioned study, the assumption of constant IGM den¬ 
sity should be released and a realistic density distribution adopted. In this case, 
the main complication would be to properly follow the propagation of ionizing 
photons into the IGM. To this aim, Ciardi et ah (2001) have developed the code 
CRASH, a Monte Carlo approach to the radiative transfer of ionizing photons 
produced by a source with known emission properties through a given density 
field. 

Using CRASH, we simulate the reionization process in a box of 30 Mpc co¬ 
moving (Stohr et ah; Ciardi, Stohr & White in prep.). This choice of the box 
dimension allows us to study a region large enough to be representative of the 
universe and at the same time to resolve objects with masses r\j 10® Mq, the 
ones responsible for the main ionizing photon production (see Fig. 1). The box 
is obtained from a combination of a high-resolution ACDM N-body simulation 
(Yoshida, Sheth & Diaferio 2001) and a semi-analitycal model of galaxy for¬ 
mation (Springel et ah 2001). In Fig. 2 the ionized regions produced by the 
simulation are shown, at z ~ 13 (left panel) and z ~ 10 (right panel), when the 
volume averaged ionization fraction is 0.1% and 4% respectively. 

The reionization process critically depends on the choice of the source emis¬ 
sion properties. Here, we have adopted a Larson IMF and a spectrum typical 
of PopIII stars; these assumptions give an upper limit to the ionizing photon 
production. More critical is the choice of the escape fraction of ionizing pho¬ 
tons, fesc for which we have adopted the indicative value of 20%. As fesc is a 
highly uncertain parameter, Ciardi, Bianchi &: Ferrara (2001) have tried to con¬ 
strain it for a Milky Way type galaxy via 3D numerical simulations, using the 
code CRASH to follow the photon propagation. They have compared results for a 
smooth Gaussian density distribution and an inhomogeneous, fractal one, with 
realistic assumptions for the ionizing stellar sources based on available data in 
the solar neighborhood. In Fig. 3, the evolution of fesc as a function of the total 
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Figure 3. Evolution of the escape fraction as a function of the total 
ionization rate for a Gaussian (filled triangles) and a fractal (open 
circles) density distribution. 


ionization rate, A/ly, is shown in the case of a Gaussian (filled triangles) and a 
fractal (open circles) density distribution. Values in the range 2% < /esc < 50% 
are found depending on A/ly and the density field, confirming the high uncertainty 

of fesc- 

4. Conclusions 

To properly describe the IGM reionization process, the following ingredients are 
needed: i) a reliable model for structure formation, including a self-consistent 
treatment of feedback effects and the best possible balance between large volume 
and high resolution of the simulation; ii) an accurate treatment of the radiative 
transfer of ionizing photons; hi) a better constraint on the parameters involved 
in the calculation, such as the source emission properties and the escape fraction. 
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